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Abstract 

Both quantum and classical behavior of single atomic spins on surfaces is determined by the 
local anisotropy of adatoms and their coupling to the immediate electronic environment. Yet 
adatoms seldom reside on surfaces alone and it is generally acknowledged that substrated-mediated 
interactions can couple single spins among each other impacting their magnetic behavior. Here we 
show that also magnetic anisotropy, which is usually considered to be a constant determined by the 
local crystal field, can be extremely sensitive to such interactions. By the example of Co dimers on 
Cu(001) and Pt(001) surfaces we highlight the intricate interplay of exchange coupling and magnetic 
anisotropy providing a much sought possibility to tune the latter through deliberate adjustment 
of the adatoms’ separation. As a technologically relevant implication we demonstrate the impact 
of such emergent non-local anisotropy on the hysterectic properties of single-atom magnetization 
curves. 
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For a while now the interest of surface science community has been keenly focused on 
nano-magnets 1 being of significance from both technological 2 ® and basic magnetism 5 7 
standpoints. With scanning probe microscopy having reached unprecedented ubiquity as 
a tool of choice for state-of-the-art basic surface science studies, the held of view of most 
studies has narrowed down to single adatoms or single-molecule-magnets (SMM). The 
key parameters defining both quantum and classical behavior of such systems, i.e. mag¬ 
netic moment, anisotropy energy (-Ema) and coupling to environment, can be determined 
for individual spins based on inelastic tunneling spectroscopy (IETS) 8 11 and single-atom 

magnetometrjff^HIll 

To extract the above parameters by magnetization or IETS measurements and to model 
the systems one usually resorts to simplified semi-empiric Hamiltonians 8 9 which in 
most cases do a fairly good job of describing the experimental observation. To account for 
dynamic effects (e.g. Kondo and quantum spin tunneling) or anisotropic interactions more 
elaborate formulations 15 18 have been resorted to. Yet to date, all of them have one thing in 
common - they treat magnetic anisotropy as an intrinsic constant parameter of the adatom 
or SMM, usually defined by the local adsorption geometry, the resulting ligand held and the 
spin-orbital coupling in the systeuP^®M. 

At the same time, in real life adatoms seldom reside on the surface alone. For example, 
it is well known that interatomic exchange interaction with neighboring impurities, for one, 
can strongly affect the magnetic order in the system across distances of several nanome¬ 
ters via direct!- 14 and substrate-mediated Ruderman-Kittel-Kasuya-Yosida (RKKY) 13 20 21 ^or 
super-exchange 22 24 type interactions and can be further directionally anisotropic 12 ! Re¬ 
cently, on the continuous quest for finding ways to deliberately tune magnetic properties of 
surface spins 25 ^3, it has been found that magnetic anisotropy can be modified by changing 
the coupling of the adsorbate spin to the substrate electron batl i 28|29 ( scanning probe-tip 
interaction 30 and hydrogenatiorP®H, but also by the introduction of neighboring adatoms 
in the immediate vicinity 32 via the changes in the surface relaxation induced thereby. 

Careful analysis of experimental data, however, reveals scattered hints that also interme¬ 
diate and long-range interaction of adatoms can have an impact on the magnetic anisotropies 
in the system. For example, the variation of anisotropies of Co adatoms on semi-insulating 
CuN islands on Cu(001) surfacd^ can only partially be explained by the variation of the 
coupling to the substrate across the island if the presence of other Co adatoms on the same 
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island is neglected. Also, in general, the scatter of anisotropy values tends to be higher in ex¬ 
periments with higher density of adsorbate^ 28 31 and on substrate exhibiting stronger RKKY 
properties (metallic substrates versus semi-insulating ones)^® 31 33 . More isolated impurities 
on decoupling substrate ^ 11 * 29 *, on the contrary, usually exhibit much smaller error-bars of 
anisotropy values. 

In the present letter we address this issue and show that electronic interaction between 
impurities does not only manifest itself in exchange coupling variation across distances of sev¬ 
eral nanometers but can equally strongly affect the anisotropy of individual adatoms. Using 
the prototypical system of Co adatoms on (001) surfaces of Cu and Pt and a combination of 
first-principles and Kinetic-Monte-Carlo (KMC) approaches we demonstrate how adjusting 
the distance between individual spins on the surface one can tune their local anisotropies and 
thus their response to external stimuli, the latter being illustrated by interatomic-distance 
dependent single-atom magnetization curves. 

Our calculations are based on the projector augmented wave (PAW) method as imple¬ 
mented in the Vienna Ab-initio Simulation Package (VASP) cod d 34 * 35 ! To exchange interac¬ 
tion is treated within the generalized gradient approximation in the Perdew-Becke-Ernzerhof 
form 36 . The substrate is modeled by a large (10 x 8 atoms) cubic super-cell slab having five- 
layers thickness and 16 A of vacuum space. The considered Co dimer configurations are 
shown in the inset of Figure |] A 3 x 4 A;-mesh is employed to sample the Brillouin zone 
with a plane-wave energy grid cut-off at 450 eV. First, scalar-relativistic calculations are 
carried out for the structural optimization and then fully-relativistic ones to determine the 
Ema- Spin-orbit coupling— and dipolar correctiond 38 ^ are taken into account in the latter 
calculations. Geometrical relaxations are performed until the forces acting on the atoms 
become less than 10” 2 eV/A and the change in the total energy between two successive 
electronic steps is smaller than 10“ 7 eV. The Ema of the dimers is calculated from the differ¬ 
ence in the sum of the one-electron Kohn-Sham band energies at a fixed potential (magnetic 
force theorem) 39 11 . To simulate the magnetization curves Kinetic-Monte-Carlo method is 
employed as described in Refs. 02] and 03] 

We begin by revisiting the well-known phenomenon of RKKY-mediated exchange inter¬ 
action of magnetic adatoms (in our case Co) on metallic surfaces. Figure [I] depicts the 
exchange interaction energy E ex = Ep m — -Fafm (defined as the energy difference between 
ferro- (FM) and antiferro-magnetically (AFM) coupled states of the dimer) as a function of 


3 










1 2 3 4 567 oo 

4-,- 1 ,-! rl ,- I I , i-(— 



separation d (A) 

Figure 1. Exchange interaction between Co adatoms deposited on Cu(001) (red circles) and Pt(OOl) 
(blue squares) surfaces as a function of adatom separation d. The rightmost point denoted oo 
corresponds to a non-interacting dimer or single adatom. Numbers on the top axis indicate the 
nearest-neighbor index corresponding to the numbers in the inset (here the color-code is: dark blue 
- Co, medium brown - first surface layer of Cu/Pt, light brown - second layer of Cu/Pt). 

the separation distance d between the Co adatoms on Cu(OOl) (red circles) and Pt(OOl) (blue 
rectangles) surfaces. In the case of Cu(OOl) surface, the results of our calculation are, as ex¬ 
pected, in good semi-quantitative agreement with the results obtained in earlier non-relaxed 
calculations based on the Korringa-Kohn-Rostoker (KKR) Green function method 20 . In¬ 
deed, structural relaxations are significant up to second-nearest neighbors (NN) (d ~ 3.7 A) 
only. For compact dimers the Co-Co bond length is significantly contracted with respect to 
the surface inter-site distance. The Co-Pt distance is also smaller than Co and Pt interlayer 
distance (about 17% relaxation for Cu and 27% for Pt). However, the in-plane relaxations 
are negligibly small for larger interatomic separations and even the relaxations towards the 
surface do not show any interplay with the observed magnetic order (which also explains 
the surprising match of our results and KKR calculations). 

A predominant FM coupling is obtained for most separations except a region of 5 < d < 
7 A where a FM-to-AFM transition is observed 20 . All in all, the exchange interaction depen¬ 
dence on the separation distance d shows a non-monotonous behavior in line with the RKKY 
physic^ 0 . For Co atoms on Pt(OOl) surface we obtain a similarly non-monotonous exchange 
interaction dependence on d (blue squares in Figure [Tj) , which is comparable in strength to 
that on Cu(OOl) at larger separations and is weaker than on Cu(OOl) for compacter dimers 
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Figure 2. Dependence of local per-atom magnetic anisotropy of (a) a Co dimer on a Cu(001) sur¬ 
face considering both ferromagnetic (FM, red circles) and antiferromagnetic (AFM, blue triangles) 
ordering of spins and the same for (b) a Co-dimer on a Pt(001) (FM in blue squares and AFM in 
orange rhombs). Positive F?ma values stand for out-of plane Ema- The rightmost point denoted oo 
corresponds to a non-interacting dimer or single adatom. The spin polarization of electrons induced 
in the Pt(001) surface by the presence of a ferromagnetic (c) and antiferromagnetic (d) Co dimers 
at 6.28 A separation. 

due to the larger Pt lattice constant. 

Conventionally, knowing the d-dependent exchange coupling strengths presented above 
along with on-site anisotropy values for a single Co adatom on Cu or Pt surface would be 
enough to fully describe the magnetic behavior of the system. However, if we take a look at 
the values of E M a self-consistently calculated for FM coupled Co dimers on Cu(001) surface 
presented in Figure [2|a) (red circles) we shall see that the values are far from being constant 
but rather exhibit a strong dependence on the interatomic separation d. Ema for a Co dimer 
Ema = E z — F/bond (defined as the energy difference between the dimers with spins aligned 
along the normal to the surface and the bond axis) is non-monotonously changing in the 
range of 0.1 — 0.6 meV with increasing d. The only exception is the compact dimer with 
stronger Ema = 1-94 meV). 

Doubtlessly, Ema of the Co adatoms is related with local environment of the adsorbates. 
However, it proved rather difficult to put a handle on the electronic origin of Ema- The 
local density of states (LDOS) of Co adatoms exhibits slightly variations with regard to 
the magnetic order of the dimers. Only the orbital-resolved LDOS analysis is meticulous 
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enough to provide insights into both the magnetic order and distance-dependence of Ema- 
A more detailed discussion of the subject can be found in the supplement. Here we shall 
only mention that using the second-order perturbation formula®, the Ema can be split 
into out-of plane and in-plane contributions by the symmetry of angular momenta matrix 
elements between the unoccupied and occupied d-states near the Fermi level. Essentially, 
the substrate-mediated coupling of Co adatoms to each other ever so slightly alters the 
occupancies of d orbitals of the adatoms resulting in the change of local magnetic anisotropy 
magnitude and even its direction. 

This being said, it should not be any more surprising that the exchange interaction be¬ 
tween Co adatoms has a pronounced effect on their magnetic anisotropy. Indeed, within 
a generalized Heisenberg Hamiltonian model the Ema can be split into on-site and inter¬ 
site contributions and renormalized to an effective second-order magnetic anisotropy which 
contains both contributions that appear from the anisotropy of the inter-site exchange in¬ 
teractions underlining their relationship 45 16 . From the electronic point of view, on can also 
look at the problem in the frame of a simple Anderson model picture 47 . According to it, 
FM coupling results in the atomic levels of the adatoms being split in both spin channels, 
while AFM coupled impurities exhibits single-atom-like electronic structures. This leads the 
above mentioned slight changes in the occupancies and therethrough to the apperance of the 
non-local anisotropy contribution. This dependence of Ema on the exchange order is well 
traceable in our calculations. The blue triangles in Figure [2j^a) denote the d-dependent Ema 
values for AFM-coupled Co dimers on Cu(001). The trends of both FM and AFM curves 
are very close, yet the AFM curve unmistakeable exhibits less variation in Ema values in 
line with the above electronic-structure arguments. Generally, the presence of interatomic 
coupling tends to increase the Ema values per atom as compared to an isolated atom (for 
Co Ema values of 0.45 meV per atom in a dimer in average as compared to 0.26 meV for a 
single adatom). 

The dependence of anisotropy of coupled Co-atoms on Pt(001) surface on the separa¬ 
tion d and their magnetic ordering is even more feature-full (Figure [2](b)) . Not only does 
the anisotropy exhibit a highly non-monotonous dependence on d featuring an oscillatory 
switching between in-plane and out-of-plane anisotropies. It also shows a much stronger 
dependence on the magnetic ordering of the adatom spins, a change of coupling sign usually 
bringing about the change of the sign of Ema- These contrasting results reveal the different 
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underlying adatom-substrate interaction mechanisms. Essentially, the interplay between 
surface-mediated and strong spin-orbit interactions (existing in the later case) together with 
the magnetic order of the dimers determine the Ema nature. First and foremost, the strong 
hybridization between the Co-3 d with the Pt-5d states results in the larger values of Ema 
for Co/Pt in comparison with Co/Cu. The striking susceptibility of anisotropy of Co dimers 
on Pt(001) as compared to Cu(001) can be understood if one considers that Pt-surface is a 
much more polarizable substrate of the two, as evidenced, f.e., by the strong directionality 
of RKKY thereon 12 and long distance character of the interaction even in absence of, e.g., a 
surface stat^- 1348 . The difference in dimer-induced substrate polarizations is shown in Fig¬ 
ure § c) and (d) for FM and AFM coupled Co dimers at d = 6.28 A. While on Pt(001) Co 
atoms show show robust and large magnetic moments of /j,qo = 2.2 /i# and a spreading sub¬ 
strate polarization cloud, on Cu(001) the magnetic moment of Co adatoms is fiQo = 1-9 hb 
and no induced polarization in the surface was observed. 

It is thus clear that anisotropy is not, as it is often believed, an intrinsic constant of 
a single adsorbate and should be treated with care. While it cannot be unambiguously 
claimed, that interaction with surrounding magnetic adatoms is solely responsible for the 
dispersion of experimentally estimated anisotropy valuers, it should definitely be duly taken 
into account when interpreting the results of past and future experiments. 

Another thing worth mentioning is that many experimentally observable magnetic char¬ 
acteristics of single adatoms and atomic assemblies 3 14 are intimately linked to such param¬ 
eters as exchange interaction and anisotropy. To illustrate the implication of the dependence 
of magnetic anisotropy in surface structures on the interatomic interactions we simulate the 
response of Co dimers at different interatomic separations d to the applied external magnetic 
held B as described by single-atom magnetization curves 13 —. The simulation is carried out 
using the stochastic Kinetic-Monte-Carlo approach which is known to yield accurate re¬ 
sults in good agreement with experimental observations 4213 . The basic characteristics of 
the system hereby are the remanence and coercivity. They are also tightly correlated with 
relaxation times of spin systemdP such as single adatoms and SMMs. 

To keep our simulations as close as possible to those used for fitting experimentally 
observed magnetization curves we describe our Co dimers on Cu or Pt surfaces with a 
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Figure 3. Hysteresis loops ( M(B)/Mg ) for representative Co dimers on Cu (a) and Pt (b) surfaces. 
The magnetization M is normalized by its saturation value Mg. 


classical Heisenberg Hamiltonian 


H = J(d ) S 1 S 2 - Y, K i( d ) S L -d Y, BSi, 

i=l,2 i=l,2 


( 1 ) 


where the exchange coupling J(d) and Ema constant K(d ) are extracted from our ab intio 
calculations. For simplicity, the anisotropy is considered to be uniaxial and the external 
magnetic field B to be always collinear to the easy anisotropy axis. The magnitude of the 


spins 


Si 


is considered constant with respect to interatomic separation and correspond to 
a magnetic moment of mco = 1-92 /xb as yielded by our first-principle calculations. The 
temperature of the simulation is taken to be 0.4 K which is within the typical experimental 
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Figure [3] shows magnetization curves for several representative Co dimers on Cu(001) 
(a) and Pt(001) (b) substrates. The magnetic field is swept from — B 0 to +B 0 and back 
(Bq = 2T(10T) for Co dimers on Cu(Pt), respectively) by increments of 1 mT making full 
cycles at a sweeping rate of dB/dt = 130 T x sec^ 1 (according to Ref. 03j). On Cu(001) sur¬ 
face Co dimers exhibit hysteresis loops with relatively small coercive fields of 87 and 120 mT 
and remanences of about 0.67 /xb for d = 3.41 and 8.11 A, respectively. It is apparent, how- 














ever, that the hysteresis walls are differently sloped hinting at a variation of the defining 
magnetic parameters. The dimer with d = 5.17 A, on the contrary, displays a paramagnetic 
response. This behavior is consistent with the idea of the first and second dimers having a 
combination of stronger exchange coupling and weaker per-atom anisotropies, or vice-versa, 
strong anisotropy but weak exchange coupling (see Figures [I] and |2](a)) , which defines the 
shapes of the hysteresis loops. The third dimer with negligible exchange, though relatively 
strong anisotropy, behaves as two independent adatoms, which are known to exhibit param¬ 
agnetic behavior due to quantum tunneling if magnetization 42 . A similar behavior is found 
in the case of the Co dimers on Pt. If the exchange interaction is small (~ 3 meV), no mag¬ 
netic response is found as in the case of the dimer with d = 5.63 A (Figure [3](b) ). Moreover, 
one observes that the larger Ema values obtained for these systems prompt higher coercive 
fields (up to 4.4 T). One can again conclude then that strong exchange interactions are re¬ 
sponsible for the hysteresis stabilization while Ema controls the shape of the magnetization 
curves. The interplay of Ema and exchange interaction defining the shape of the hysteresis 
loops of single atom magnetization curves can be of import when analyzing the the results 
of corresponding experimental studied®, though in many cases the impact of the inter¬ 
atomic interaction shall be somewhat weaker, especially when the adatoms are placed on 
semi-insulating substrates or more sparsely distributed on a surfaceP^ES. 

To draw a bottomline, with our study we underline the importance of departing from the 
oversimplified picture of anisotropy being a purely local constant and considering the ap¬ 
parent influence thereon of the interatomic interaction present especially in ensembles of 
adatoms deposited on substrates exhibiting RKKY features. This generally non-trivial sep¬ 
aration dependence of the anisotropy can have a appreciable effect on such characteristics of 
magnetic systems as Kondo behavior, quantum spin tunneling and the resulting observable 
quantities as the inelastic tunneling spectra and the single-atom magnetization curves. 
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